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We have developed a new method of an earthquake- resist ant design to support conventional 
aseismic designs using acoustic metamaterials. We suggest a simple and practical method to reduce 
the amplitude of a seismic wave exponentially. Our device is an attenuator of a seismic wave. 
Constructing a cylindrical shell- type waveguide that creates a stop-band for the seismic wave, we 
convert the wave into an evanescent wave for some frequency range without touching the building 
we want to protect. 
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I. INTRODUCTION 

Earthquakes are the result of sudden release of huge 
amount of energy in the Earth's crust that produces 
seismic waves. A sudden forthcoming of the seismic 
waves with large amplitudes and low frequencies have 
been great hazards to life and property. It is the col- 
lapse of bridges, dams, power plants, and other struc- 
tures that causes extensive damage and loss of life during 
earthquakes. Aseismic capabilities are highly relevant to 
public safety and a large amount of research has gone 
into establishing practical analysis and design methods 
for them. Numerous earthquakeproof engineering meth- 
ods have been tried to resist earthquakes, but still we are 
not so safe. 

Seismic waves are a kind of inhomogeneous acous- 
tic wave with various wavelengths. There are two 
types of seismic waves: body waves and surface waves. 
P (Primary) and S (Secondary) waves are the body waves 
and R(Rayleigh) and L(Love) waves are the surface 
waves. Surface waves travel slower than body waves and 
the amplitudes decrease exponentially with the depth. It 
travels about 1 ~ 3km/ sec with lots of variety within the 
depth of a wavelength Q . The wavelengths are in the 
order of 100m and the frequencies are about 10 ~ 30Hz^ 
that is, low end and just below the audible frequency. 
However, they decay slower than body waves and are 
most destructive because of their low frequency, long du- 
ration, and large amplitude. 

Rayleigh waves can exist only in an homogeneous 
medium with a boundary and have transverse motion, [ll, 
Q Earthquake motions observed at the ground surface 
are mainly due to R waves. On the other hand, L waves 
are polarized shear waves guided by an elastic layer. It 
is this that causes horizontal shifting of the Earth during 
earthquakes. L waves have both longitudinal and trans- 
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verse motion and this is what most people feel directly 
during earthquakes. 

Recent development of metamaterial science opens a 
new direction to control the seismic waves. Farhat et al. 
proposed a design of a cloak to control bending waves 
propagated in isotropic heterogeneous thin plates [s'-dj. 
Their cloaking of shear elastic waves pass smoothly into 
the material rather than reflecting or scattering at the 
material's surface. However, the cloaked seismic waves 
are still destructive to the buildings behind the cloaked 
region. 

In this paper we introduce a method to control seismic 
waves by using a new class of materials called metamate- 
rials. Metamaterials are artificially engineered materials 
which has a special property of negative refractive index. 
We present here a solution that a metamaterial acts as 
an attenuator by converting the destructive seismic wave 
into an evanescent wave by making use of the imaginary 
velocity of stop-band of the wave. 

There are many representations to express the scale 
of earthquakes. Among them the magnitude that comes 
from the amplitude of the seismic waves is most impor- 
tant. The common form of the magnitude, M, is the 
Richter-scale defined by comparing the two amplitudes 
in a logarithmic scale as 



M = log— , (1) 



where A is the maximum amplitude of the seismic wave 
and Ao is the maximum amplitude of the background 
vibration and order of /im. The equipment measures a 
transformed magnitude of the intensity. Another factor 
of the strength is PGA(Peak Ground Acceleration) but it 
is not expressed in a closed form. We focus here on how 
to reduce the amplitude of the seismic wave by using the 
properties of the metamaterials. 
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II. NEGATIVE MODULUS 

Acoustic waves are created by compressibility or elas- 
ticity of the medium. Young's modulus, F, is a one- 
dimensional compressibility defined by AP = YAl/l^ 
where AP is the pressure or stress and / is the length. 
Shear modulus, G, is a two-dimensional one for a surface 
wave defined by AP = GAx/h, where Ax is the horizon- 
tal shift and h is the height of the object. Bulk modulus, 
P, is a three-dimensional one for a body wave defined by 
AP = -BAV/V. 

Seismic medium of Earth crust can be considered as 
an accumulation of infinite number of elastic plates. Al- 
though the seismic surface wave is not pure two dimen- 
sional, the velocity is mainly dependent upon the density, 
p, and shear modulus, G, of the seismic medium. Seis- 
mic wave is a kind of acoustic wave and every acoustic 
wave propagates following two wave equations in princi- 
ple. Assuming the plane wave time dependence e*^^ the 
pressure, p, and the velocity, v of the wave in the two 
dimensions are expressed as the Newton's 2nd law 



and the continuity equation 

iujp — GVs • iT, 



(2) 



(3) 



where Vg is the Laplacian operator at the surface, p is 
the pressure, uj is the anguar frequency of the wave, and 
V is the velocity. 

The Eq. ([2]) and Eq. (|3|) generates the wave equation 

as 



V^p + — p = 0, 
where the velocity of the seismic wave is 



(4) 
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If the shear modulus becomes negative, the velocity be- 
comes imaginary. Then, so does the refractive index, n, 
or the inverse of the velocity as n = Vq/v = Vo^/p/G, 
where Vq is the background velocity. Since k = 27rn/A, 
the wave vector becomes imaginary, too. Therefore, the 
imaginary wavevector makes the amplitude of the seismic 
wave become an evanescent wave. We call it a barrier or 
attenuator. Note that the impedance Z = pv = ^fpG 
becomes imaginary because it is an absorption. 

Negative shear modulus of elastic media has been stud- 
ied and realized very recently [6|-[8|. The key was the 
Helmholtz resonators similarly with the negative bulk 
modulus. The Resonance of accumulated waves in the 
Helmholtz resonator reacts against the applied pressure 
at some specific frequency ranges. Then, the nega- 
tive modulus is realized by passing the acoustic wave 
through an array of Helmholtz resonators. Therefore, 
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FIG. 1: Real and imaginary parts of the effective shear modu- 
lus. The negative peak of the imaginary part means an energy 
absorption. 7 = F/cJo- 



the sound intensity decays exponentially at some reso- 
nant frequency ranges. In general the elastic material 
is described by three independent effective parameters of 
G, P, and p. Therefore, sometimes the G is replaced bya 
linear combination of G, P, and the Lame constant [l|,[2|], 
but it will not change the structure of the theory. Acous- 
tic waves from the modulus share fundamental properties 
of sound waves. 

From the formalism of electromagnetic response in 
metamaterials, effective electric permittivity and effec- 
tive magnetic permeability show negative values at some 
specific frequency ranges around resonances [9]. The 
Helmholtz resonator is a realization of an electrical reso- 
nance circuit by mechanical correspondence. It is known 
as that the plasmon frequency in metals or in an array 
of metal wires produces the electric permittivity as [HI 



uj(uj + zF) 



(6) 



where cjp is the plasma frequency and F is a loss by damp- 
ing. 

The Eq. (j2j) is the counterpart of the Faraday's law and 
Eq. (|3]) is of the Ampere's law by the analogy of electro- 
magnetism and mechanics. The inverse of the modulus 
in mechanical system corresponds to the electric permit- 
tivity in electromagnetic system. Considering the struc- 
tural loss, the general form of the effective shear modulus, 
Ge/ / , is given similarly with the general form of the bulk 
modulus as 
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where is the resonance frequency and P is a geometric 
factor p^, ll6| . The real and imaginary part of the Ge/ / 
is plotted in Fig. [T]at specific values of P and F. The real 
part can be negative at resonance and slightly increased 
frequency ranges. The negative range of the real part is 
the stop-band of the wave. When the imaginary part, the 



loss, is small compared with real part, the effective shear 
modulus has negative value dX 1 < uo /uOq < Vl + F. 



III. SEISMIC ATTENUATOR 

We can built an attenuator or an earthquakeproof bar- 
rier of a seismic wave by filling-up many resonators un- 
der the ground around the building that we want to pro- 
tect. Then, the amplitude of the seismic wave that passed 
the waveguide is reduced exponentially by the imaginary 
wavevector at the frequency ranges of negative modu- 
lus. Mixing up many different kinds of resonators will 
cover many different corresponding frequency ranges of 
the seismic waves. 

If we assume that the plain seismic wave of wavelength 
A propagates in x— direction, the amplitude of the wave 
reduces exponentially as 



Ae 



ikx 



(8) 



Let the initial seismic wave, that is, before entering the 
waveguide, have amplitude Ai and magnitude M^, and 
final seismic wave, that is after leaving the waveguide, 
have amplitude Af and magnitude Mf following the Eq. 
([1]). Then, Af is written as Ai from Eq. ([8|) as 



^.g-27r|n|x/A ^^^^ 



(9) 



The amplitude of the seismic wave reduces exponen- 
tially as passing the waveguide of metamaterials. We 
can rewrite Eq. (|9|) with the definition of the magnitude 
in Eq. ^ as 



(10) 



Taking logarithms both sides of Eq. (p^Oj) . we obtain the 
width of the waveguide, x Ax, as 

Ax —— = , , AM, (11) 

27T \n\ \n\ ' ^ ^ 

where AM = Mi- Mf. 

For example, if the refractive index is n = 2 and the 
wavelength of the surface wave is A = 100m, we need the 
waveguide of the width Ax 18m to reduce AM = 1. 
If the aseismic level of the building is M = 5 and the 
width of the waveguide surrounding the building is about 
60m, then the effective aseismic level of the building is 
increased to M = 8. Therefore, a high refractive index 
material is desirable for a narrow waveguide. In civil 
engineering earthquakeproofing methods must be practi- 
cal, that is, clear to manufacture and easy to construct. 
The resonator should be easy to build. We designed an 
example of a resonator in Fig. [2l 

The size of the cylinder can be estimated from the 
analogy between electric circuits and mechanical pipes. 
A pipe or tube with open ends corresponds to an induc- 
tor, and a closed end corresponds to a capacitor [H, [I^- 
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FIG. 2: (a)A sample of a meta-cylinder with 4 side holes. 
The size of the cylinder is less than the wavelength of the 
surface waves, (b) A combined form of the 4 meta-cylinders. 
An electrical analogy is shown. 



where p is the density inside the volume, V is the effective 
length, S is the area of the cross-section, V is the volume, 
and V is the velocity inside. From Eq. ([12]) the resonant 
frequency is given as 




(12) 



(13) 



In the meta-cylinder V is the effective length which is 
given by / + 0.85(i [17], where / is the length of the 
hole or thickness of the cylinder, and d is the diameter of 
the hole. 

An example of the design of the meta-cylinder for the 
seismic frequency range is followings: the diameter of the 
hole is order of 0.1m, the thickness of the cylinder is or- 
der of 0.1m, and the volume inside is order of m^. Since 
the meta-cylinder is considerably smaller than the corre- 
sponding wavelengthes, the array meta-cylinders behaves 
as a homogenized medium. 

The shape of the meta-cylinder is neither necessary to 
be circular nor to have 4 holes. It could be any form of 
a concrete box with several side holes. Cubic or hexago- 
nal boxes would be fine. Various kinds of resonators may 
cover various kinds of resonance frequencies of the seis- 
mic waves. There happens an energy dissipation of the 
seismic waves inside of the waveguide and the absorbed 
energy will turn into sound and heat. They makes the 
temperature of the waveguide increasing depending on 
the magnitude of energy that arrives at the waveguide. 

A vertical view of the metamaterial barrier with many 
meta-cylinders are in Fig. [3l The width of the barrier. 
Ax, is predicted in Eq. (fTTj) . The depth of the waveguide 
should be at least the foundation work of the building to 
protect as in Fig. [3l but it is not necessary to be more 
than the wavelength of the surface waves. The completed 
form of the waveguide is the aseismic cylindrical shell of 
many concentric rings in Fig. HJ 

Seismic waves cannot pass through water. Then, we 
can imagine a water barrier of a big trench filled with wa- 
ter. However, it is not stable to stand a series of attacks 
of 'foreshock main shock aftershock.' Because, if 
the outer part of the water trench is brought down by a 
seismic wave, we do not have enough time to rebuild it 
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FIG. 3: A vertical landscape of the metamaterial(MTM) 
barrier and the building to protect from the seismic wave. 




FIG. 4: A sky view of a completed seismic waveguide with 
many meta-cylinders. 



to prepare for the main shock and the after shock. Wa- 
ter cannot have a high refractive index n in Eq. (pT]) , 
too. Maintaining the depth of the water trench up to the 
wavelength of the seismic wave require huge amount of 
water and, therefore, not be practical. 



IV. SUMMARY 

We introduced a supportive method for aseismic de- 
sign. It is not to add another aseismic system to a build- 
ing but to construct an earthquakeproof barrier around 
the building to be protected. This barrier is a kind of 
waveguide that reduces exponentially the amplitude of 
the dangerous seismic waves. 

Controlling the width and refractive index of the 
waveguide, we can upgrade the aseismic range of the 
building as needed in order to defend it, at will, without 
touching it. It could be a big advantage of the waveg- 
uide method. This method will be effective for isolated 
buildings because we need some areas to construct the 
aseismic shell. It may be applicable for social overhead 
capitals such as power plants, dams, airports, nuclear re- 
actors, oil refining complexes, long-span bridges, express 
rail-roads, etc. 
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